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1. Introduction and Motivation 

In many condensed matter systems, one finds phase transitions governed by fixed 
points which exhibit "dynamical scaling" 

t -> X z t, x -> Ax, z^l (1.1) 

instead of the more familiar scale invariance which arises in the conformal group 

t -»• At, x -> Ax . (1.2) 

A toy model which exhibits this scale invariance (and which is analogous in many ways 
to the free scalar field example of a standard conformal field theory) is the Lifshitz field 
theory: 

C = J d 2 xdt {{d t (j)f - k(V 2 <P) 2 ) . (1.3) 

This theory has a line of fixed points parametrized by k and arises at finite temperature 
multicritical points in the phase diagrams of known materials [QJJ] . It enjoys the anisotropic 
scale invariance ( |1 . 1| ) with z = 2. 

This fixed point and its interacting cousins have become a subject of renewed interest 
in the context of strongly correlated electron systems. For instance, in the Rokhsar- 
Kivelson dimer model 0, there is a zero-temperature quantum critical point which lies 
in the universality class of ( |1.3| ) |4|] (for a nice general exposition of the importance of 
quantum critical points, see ||). Similar critical points also arise in more general lattice 
models of strongly correlated electrons fj,[jj§]. The correlation functions in these models 
have interesting properties like finite-temperature ultra-locality in space at fixed time 0, 



which may be important in explaining certain experimental results [ ITD ,|TT ] . Such theories 
are also of interest in 1+1 dimensional systems ||T2||T3"| . 

Furthermore, such fixed points seem to have a non-trivial generalization to non- 



Abelian gauge theories. The Lagrangian (|1.3|) can be dualized to that for an Abelian 



gauge field in a standard way, since scalars are dual to vectors in 2+1 dimensions. This 
yields a Lagrangian with the unusual property that the usual E 2 term has vanishing co- 
efficient; the leading terms in the Lagrangian depend on derivatives of the electric field. 
Freedman, Nayak and Shtengel analyzed a similar theory with SU(2) gauge group in |TJ 



They provided evidence that the 577(2) theory as well, has an interacting line of fixed 
points with z = 2. However, these fixed points are strongly coupled; their existence and 
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their detailed properties (correlation functions, etc.) lie outside of the regime where the 
analysis in |TJJ] was performed most reliably. 

AdS/CFT duality |15| , ^6| , |17|| provides a well known technique to obtain weakly cou- 
pled and calculable dual descriptions of strongly-coupled conformal theories, in terms of 
gravity or string theory on a weakly curved spacetime. The symmetries of the gravitational 
background geometrically realize the symmetries of the dual field theory; so for instance 
the conformal group SO(D, 2) of a D-dimensional CFT, arises as the group of isometries of 
AcISd+i- It is natural to ask: can we find a more general class of spacetimes which could 
be dual to theories with non-trivial dynamical critical exponents as in (|1 . 1|) ? The general 
such theory has far fewer symmetries than a conformal theory; it enjoys scale invariance, 
invariance under spatial and temporal translations, spatial rotation invariance, and P and 
T symmetry.0 

In this paper, we attack this question. While we will focus most concretely on the 
case z = 2 (where many of the most interesting examples arise in the condensed matter 
literature), the techniques we use are clearly more general. In Sec. 2, we show that 4D 
gravity with a negative cosmological term, in the presence of a modest set of p-form gauge 
fields, can support metrics which geometrize the symmetries of these Lifshitz-like fixed 
points. In Sec. 3, we compute two-point correlation functions for the simplest scaling 
operators, those dual to free bulk scalar fields. Unlike the situation in conformal field 
theories, these correlators contain more information than can be inferred from the scaling 
dimension of the operator alone. We study the renormalization group flows between our 
fixed points and conventional conformal field theories (analogous to the flow that would 
arise by perturbing the Lagrangian ( |1.3| ) by the operator — (V</>) 2 ) in Sec. 4. Finally, we 
conclude with a discussion of several interesting questions that we hope to address in the 
future. 



1 Non-relativistic theories without particle production and with dynamical scaling enjoy a 
larger symmetry group. When z = 2, the enhanced symmetry group is called the Schrodinger 



group. Gravity duals for such theories have been studied in the recent works [18], initiated by the 
papers of Son and of Balasubramanian and McGreevy. One potential application seems to be to 
systems of cold atoms at the unitarity limit. The theories we study, while also lacking Lorentz 
invariance, have particle production; their structure and their potential applications are quite 
different from those in [18|. Other recent papers applying AdS/CFT duality to study different 
problems of condensed matter physics in a similar spirit have appeared in e.g. [1£]. 
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2. Gravitational Solutions 

In this section, we find gravity solutions which have the right properties to be dual 
to (interacting generalizations of) Lifshitz fixed points. We work with general z here, but 
specialize to z = 2 in the next two sections when we compute correlation functions and 
study renormalization group flows. 

2.1. The Metric 

We would like a metric invariant under the modified scale transformation ( |1 . 1| ) ■ In 
addition, we wish to study field theories invariant under time and space translations, spatial 
rotations, spatial parity, and time- reversal. We assume, as in AdS/CFT, that the "scale" 
in the dual field theory is geometrized by the presence of an additional radial dimension 
on the gravity side of the duality, and that rescaling of this radial coordinate geometrizes 
the scale transformations (|1 . 1| ) . 

These assumptions lead us to the family of geometries (one for each value of z): 

ds 2 = L 2 (-r 2z dt 2 + r 2 rfx 2 + ^ , (2.1) 

where < r < oo, <ix 2 = dx\ + ... + dx\, and L sets the scale for the radius of curvature 
of the geometry. We set the Planck length l p \ = 1, and hence every quantity above is 
dimensionless. The scale transformation acts as 

t \ z t, x -> Ax, r -> y • (2-2) 

A 

z = 1 gives the usual metric on AdSd+2, which enjoys the larger symmetry SO(d +1,2). 

This metric is nonsingular. All local invariants constructed from the Riemann tensor 
remain finite everywhere, and in fact are constant. The latter fact follows from the sym- 
metries of the geometry. The Ricci scalar takes the value —2(z 2 + 2z + 3)/L 2 . The metric 
is not geodesically complete and has peculiar behavior near r = 0; for z ^ 1, although the 
curvature invariants are small for sufficiently large L, an infalling object experiences very 
large tidal forces. This behavior is familiar from the solutions corresponding to various 



string theoretic brane systems ||20|| . The lesson from these systems is that a metric like 
( | . 1|) is physically sensible if there exists a regular black hole solution that approaches it 
in an extremal limit; we leave the construction of such black holes to future work.i 

Henceforth, we will focus on the case d = 2 (appropriate for gravitational duals to 
2+1 dimensional field theories); so we will be studying gravity in four dimensions. 



We are grateful to Sean Hartnoll and Gary Horowitz for bringing this subtlety to our 
attention. 
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2.2. Full Solution 



We would like to obtain the metrics (|2.1|) as solutions of the field equations of General 
Relativity coupled to some matter content (which can presumably arise in the low-energy 
limit of string or M-theory). Because our theory, like the toy model ( |1.3| ), may be expected 
to flow to normal CFTs under relevant perturbations, a good starting point will be gravity 
with a negative cosmological term; this will be able to capture the end of any such RG 
flows via AdS/CFT duality. 

Einstein gravity with a negative cosmological constant alone does not support the 
metrics ( j2.1| ). However, string theory also generically gives rise to p-form gauge fields. A 
modest choice of such additional content, which can support the metrics ( |2.1p , involves the 
addition of gauge fields with p = 1, 2. 

The Lagrangian is then given by: 



S= J d 4 xV=g (R - 2A) - \ f (^F{2) A*F (2) +F (3) A*F {3) ^-cJ B (2) AF {2) . 



(2.3) 



where R(2),(3) are the field strengths for the gauge fields, F( 2 ) = dA^,F^ = dB^y In 
addition to the standard Einstein-Hilbert action, a 4D cosmological constant A, and the 
kinetic terms for the gauge fields, we introduced a topological coupling between the two 
and three form fluxes f B^) A Fp) = J -4(i) A Fr$\ (up to boundary terms). This coupling 
is necessary to find appropriate solutions of the equations of motion. The topological 
coupling c needs to be quantized in many (but not all conceivable) cases, as discussed in 
detail in for instance appendix A of [21]. However, note that after redefining F( 2 ) ~~ > e^( 2 )' 



the action can be written with no e in front of the gauge kinetic term, but with c — > ce; 
in this way, we can consider c to be arbitrarily small (at weak gauge coupling), and set 
e = 1. We proceed with this convention henceforth. 

To source the metric ( |2.1| ), we need to turn on background two and three- form fluxes 
that preserve the symmetries: 

r {2) = a 9 r a e u f (3) = b e r a e x a e y . (2.4) 

Here we work with a non-coordinate basis for the one-forms: 

d,v 

t = Rr z dt, d xi = Lrdx\ 9 r = R— . (2.5) 

r 



In this basis, the metric simplifies to diag(— 1, 1, 1, 1). We fix our convention to be 
eg r g t g x g y = 1. Both fluxes are then closed. Their field equations 



d * F {2) = -cF (3) , d * F (3) = cF (2) (2.6) 

then fix the value of the dynamical exponent z in terms of the topological coupling c and 
the radius of curvature L of the dual geometry: 

2z = (cL) 2 . (2.7) 

The ratio between the two fluxes is also fixed in terms of the dynamical exponent: 

-B = Vf ' (2 ' 8) 

Note that for small c, we can obtain a weakly curved geometry with any desired value of 
z. 

The Einstein equation for the action ( |2.3| ) 

(2.9) 



G + Aa = — I vF F p2 " Pp — -a F fp^-Pp) 

p=2,3 P ' ^ ' 



then fully determines the required values of the background fluxes and the cosmological 
constant. As a result of the symmetry of ( |2.1| ) and ( ^.4| ), and the Bianchi identity, there 
are only two independent equations from the Einstein equations, which can be taken to be 
the tt and xx components: 

z 2 + 2z A 2 B 2 



L 2 ' I 2 

z + 2 A 
A + 



2 (2-10) 



L 2 4 

^) and ( |2.10[ ) then determine the necessary values of the cosmological constant and the 
fluxes: 

. z 2 + z + 4 

A = n , 

2L 2 ' 

A 2 = 2Z{Z L ; 1] , (2.11) 

2 4(z-l) 
B = L 2 ' 

We notice that reality of the fluxes requires z > 1. The dual field theories may exhibit 
"critical slowing down," but never "critical speeding up." 
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3. Two-Point Correlation Functions 

We now explore the boundary observables defined by the theory in the bulk (|2.1|) by 
generalizing the usual holographic dictionary. We specify henceforth to the case z = 2, for 
illustration. There are a variety of boundary observables one could consider. We focus on 
the simplest possibility, two point correlation functions of scalar operators, in the present 
note. After discussing generalities about the solutions for a bulk scalar field, we calculate 
the two-point function for a marginal operator (dual to a massless scalar) in detail in 
Sec. 3.2. We briefly discuss the general results for operators dual to massive scalars in 
Sec. 3.3. We note that because the symmetries in our theories are much less constraining 
than the full conformal group, the general two-point function includes an unknown scaling 
function of x 2 /|t|, and is not determined entirely by symmetry alone; so the two-point 
function already contains more non-trivial information than just a scaling dimension. 

3.1. Real scalar field 

Consider a free, real scalar field <fi in the Euclidean version of the background metric 
( p7T|) . We work in the coordinate it = 1/r so that the boundary is located at u = 



In a string construction such scalars could arise from the moduli of the compactification 
manifold, from Kaluza-Klein modes of the metric and p-form fields, or from excited string 
states. The metric fluctuation u 2 Sg xy , among other things, also satisfies the same equation 
of motion as a bulk scalar with m 2 = along the radial and temporal directions. 
The field equation for <fi is, 





The scalar has an action, 




(3.2) 



d 2 u( f>--d u( p + u 2 d 2 cp+(d 2 x + d 2 y 



)0 




= 0. 



(3.3) 



Near the boundary a solution takes the asymptotic form, 




(3.4) 
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where A±, A + > A_, are the two roots of the equation 



A(A- 4) =m 2 L 2 . (3.5) 

The requirement that the Euclidean action be finite imposes bounds on the allowed 
values of A and hence the mass of the field as in the case of the usual AdS/CFT corre- 



spondence |[22|| . For 

m 2 L 2 > -3 (3.6) 
there is a unique choice of boundary condition 

<f>(u, r, x, y) -> u A+ (0(r, x, y) + 0(u 2 )) (3.7) 

and via this choice, the scalar field is dual to an operator of dimension A + > 3. In the 
window, 

-4 < m 2 L 2 < -3 , (3.8) 

in addition to the above choice which remains valid, one can modify the Euclidean action 
to 

~ J d 4 Xv^(/>(-V 2 + m 2 )^ (3.9) 
by subtracting an infinite boundary term, so that the boundary condition 

<f>(u, r, x, y) -> u A ~ (0(r, x, y) + 0(u 2 )) (3.10) 

also leads to finite action. Thus there are two different quantizations for the scalar field 
in the range (^^), and correspondingly there are two different non-Lorentz invariant fixed 
points on the boundary: one with an operator of dimension A + > 2, the other with an 
operator of dimension 1 < A_ < 2. For m 2 L 2 < —4, the theory has a real instability; this 
is the analogue of the Breitenlohner-Freedman bound [53[ for the model in 2+1 dimensions 



and with z = 2 that we are studying. The extension of these results to general spacetime 
dimension and dynamical exponent is straightforward. 



3.2. Correlation functions of marginal scalar operators 

As in standard AdS/CFT, boundary correlators are given by the value of the renor- 
malized bulk action for specified boundary values of the bulk field (up to cut-off dependent 
field renormalizations) . To illustrate the idea, consider first a free, real scalar field with 
m = for which no field renormalization is necessary. This corresponds to an exactly 
marginal operator with scaling dimension four if we ignore gravitational loop corrections. 
(Remember that the dual field theory lives in 2 + 1 dimensions with dynamical exponent 
z = 2). 

To compute the action, we find the boundary to bulk propagator G(u, x; 0, x') 
that gives the nonsingular bulk field configuration 0(it, x) for any smooth boundary value 
</>(0, x) (x denotes collectively r, x) 

<p(u, x) = J d 3 x' </>(0, x 1 ) G(u, x; 0, x'). (3-H) 

The translational invariance in r and x makes it easy to work in the Fourier space (k 
denotes collectively us, k) where Q3.11J ) becomes 

0(u, k) = G(u, k) 4>(0, k) (3.12) 

and G(u, u, k) is the solution of (OP with m = in Fourier space 



dlG - -d u G - (u 2 z 2 + \k\ 2 )G = 0. (3.13) 

u 

with the boundary conditions 

(5(0, u, k) = 1 (3.14) 
and G(u, k) being finite as u — > oo. These uniquely determine the propagator 

G(t*, k) = e-^' 2 r(|^ + \) _ 1 _i, (3.15) 



T v 4|w| 2' 

where U(a, 6, z) is the confluent hypergeometric function of the second kind. Note that 
G(u, k) vanishes as u — * oo. 

By standard integration by parts, an on-shell bulk action is determined by the values 
of the field on the boundary 

/poo 
d 3 x J du (-^y/gg^d^ + d^y/gg^^)) 

= J d 3 x[^g uu (Pd u( p]™ (3.16) 

= J d 2 kdu(<p(0, k, u) ^(kjw) 0(0, -k, -u)) . 



We have cut off the whole space at u = e to regulate the bulk action. The 'flux factor' T 
is 

" _ _ "1 OO 

^(k,w)= G(u, -k, -u)y/gg uu d u G{u, k, u) . (3.17) 



Since the propagator G vanishes at u = oo, JF only receives a contribution from the cutoff 
at u = e. The momentum-space two-point function for the operator 0^ dual to <fi is 
calculated by differentiating ( |3.16| ) twice with respect to 0(0, k, cu), 

(<^(k, w)0^(-k, -u)) = ^(k, w). (3.18) 

Near u = 0, G has the expansion, 

G =1 - + f- [3k 4 - 20^ 2 + 2 7 (4^ 2 - k 4 ) 

4 64 L y 



+ 8k 2 M + (4^ 2 - k 4 ) log(^ 2 ) + 2(4^ 2 - k 4 )^(^ + ^-r)] + 0(u e 

2 4 w 



(3.19) 



where ij)(x) = V / (x)/V(x) is the digamma function and 7 ~ .577 is the Euler-Mascheroni 
constant. Plugging the expansion (|3.19| ) into ( |3.17| ) we pick out the leading non-polynomial 



piece in either k or w. This gives the correlation function as we take the cut-off e to 

(O^k, w)04-k, -w)) = - ^k 2 M - ^-{Au 2 - k 4 ) log \u\ 

L 2 3 k 2 (3 ' 20) 

L> , . 2 1,4 N > ' K n 



2 4w 



Specifically, the divergence arising as e — > from the term proportional to u 2 in ( |3.19| ) is 
removed via local boundary terms E3], and the terms proportional to it 4 in the first line 



of ( |3.19j ) give rise to uninteresting contact terms in spacetime. Terms 0(u 6 ) and higher 
vanish as the cutoff is removed e — » 0. 

Terms of order u 4 in the second line of (|3.19|) are the only contributors to ( |3.20|) . 
Interestingly, among the three terms in ( |3.20| ), only the last one gives rise to correlations 
between points with both spatial and temporal separation. The first two contribute terms 
localized in space. The existence of such spatially localized terms is forbidden in Lorentz 
invariant theories, and may well be related to the "ultra-local" behavior observed in H. 



Note also that (|3.20|) has the correct scaling behavior for the two-point correlator of a di- 
mension four operator. Its momentum and frequency dependence is complicated, however. 

One can understand the large distance fixed time behavior of the correlation function 
by Fourier transforming ( |3.20|) to position-space. Only the piece containing the digamma 
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function contributes in this regime, and its Fourier transform can be computed by utilizing 
the expansion 



x ' \ x + n n 

n=l 

It suffices to quote that the large distance behavior of the two point function is just a 
simple power law 

(C^(x, t)0+{% 0)> - |x| - oo. (3.22) 

This is interesting because a priori, one might have expected further suppression of the 

|x| 2 

power law result by some scaling function of x 2 /|t| (e.g. e 1*1 , which would result in an 
ultra-local equal time two-point function) could have arisen, since it is consistent with the 
symmetries of the theory. We note that the correlators of the simplest scaling operators in 
the free Lifshitz theory (|1.3|) also exhibit pure power law decay at large spatial separation 



3. 3. Correlation functions of generic scalar operators 

It is a direct generalization of the above to compute the two point functions for opera- 
tors which are dual to massive scalar fields, complicated only by additional renormalization 
effects. Specifically, to define the boundary to bulk propagator, we need to cut off the space 
and put the boundary at u = e. The propagator is now a solution to 

QlG - h u G - {u 2 z 2 + |k| 2 + ^-^)G = (3.23) 

u z l 

satisfying the boundary conditions 

G(e,k,w) = l (3.24) 

and G nonsingular for u > e. 

For generic values of m, this determines 



G(u, k) oc e -!M"VW4+OT ^ M_ + 1 + ^ + ULL 1 + v/ 4 + m 2 L 2 , \u\u 2 ) (3.25) 

with the proportionality constant set by ( |3.24|) . Now the leading non-polynomial term in 
G that contributes to the flux factor would come at "order" w 2 +v / 4+m 2 L 2 

2-V4+m 2 L 2 



G(u, k) = (^) [1 + ... + (\lo\u 2 ) VWIJ 



+ i + (3-26) 



x r(-v/4 + m 2 L 2 )r(f^ + ^ + yi + 
x r(v^TW) r ( £ + ^-AAT^) + "" 
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The two point function in momentum space is again given by the flux factor, and after 
throwing away divergent and vanishing pieces as we take e — - > 0, we find 



(C(k, u)0(-k, -w)> = -(2 + v / 4 + m 2 L 2 ) L 2 \u\ V4+m ' 2L ' 2 



F(-V4 + m 2 L 2 ) rCfe + | + ^1 + gg) (3.27) 



>4|w| 1 2 V 4 

Again ( |3.27| ) has the correct scaling behavior for the correlator of an operator of dimension 



2 + V4 + m 2 L 2 , but it also involves additional momentum and frequency dependence 



4. Holographic RG flow between z = 2 Lifshitz fixed points and z = 1 conformal 
field theories 

It is natural to expect that theories with anistropic scale invariance can flow, under 
relevant perturbation, to fixed points which enjoy full conformal invariance. For instance, 
for the toy model (|1.3|), perturbation by — (V</>) 2 induces RG flow to the free massless 
scalar field theory, which is a conformal field theory. 

Here, we find that our gravitational duals of Lifshitz-like theories come with relevant 
perturbations that induce a similar flow. We will look for solutions which represent holo- 
graphic RG flow from the z = 2 Lifshitz-like fixed point in the UV (at large r) towards an 
AdS^-like spacetime in the IR (small r). A metric ansatz which is sufficiently general to 
capture this flow is: 

ds 2 = L 2 {-r^f(r) 2 dt 2 + r 2 (dx 2 + dy 2 ) + g(r) 2 dr 2 /r 2 ) . (4.1) 

If f(r) = g(r) = 1, this is the spacetime dual to the z = 2 fixed point. If f(r) is 
inversely proportional to r, and g(r) is constant, it is AdS^. Introduce again the non- 
coordinate basis of 1-forms in which the metric simplifies to constant Minkowskian form 

9 t = Lr 2 f(r)dt 
9 X = Lrdx 
6 y = Lrdy 
9 r = Lg(r)dr/r 



Note that for non-generic values of the scalar mass m L = n — 4 with n being non-negative 



integers, the above formula ( 3.27 ) does not apply. 
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In terms of this basis, the two and three form fluxes are 

F {2 ) = ^h(r)e r A9 t (4.3) 

F {3) = ^j(r)9 r A9 x A9 y . (4.4) 

The advantage of using the 9 basis is that the fluxes remain nonsingular as long as the 
coefficient functions h(r) and j(r) are finite. In addition to the fluxes, we also have a 
negative cosmological constant A = — 5/L 2 . The z = 2 Lifshitz fixed point has h(r) = 
j(r) = 1 and the IR AdS^ has h(r) = j(r) = 

The equations of motion for the fluxes and the metric, after proper massaging, give: 

2rf'/f = (5-h 2 +j 2 )g 2 -5 
rg' = ^g*{h 2 +j 2 -h) + -g 

(4-5) 

rh! = 2gj — 2h 

rf = 2gh + -j + ^jg 2 (h 2 - j 2 - 5) . 

The counting of the equations works as follows: the two and three form fluxes each satisfy 
a single equation of motion, and the Einstein equation for our background ( fll| ) , (|4.3| ) , and 
( |4.4| ) results in only two independent equations. It takes some algebra to show the latter 
fact. It is easy to check that the z = 2 Lifshitz fixed point and AdS^ are the only two fixed 
points of the above flow equations. 

Note the relation between the above four equations. The last three form a closed set 
of ODE's for the three functions g(r),h(r) and j(r). Given any initial values for them 
at some fixed ro, we can determine their behavior as functions of r. Then by the first 
equation, we can determine the behavior of f(r). Note that the initial value of / is not 
a physical parameter, but a gauge choice. It can be changed by a global rescaling of the 
t-coordinate. 

Linearization near the fixed points gives information about the attractiveness of 
the fixed points. Around the AdS^ fixed point, the flow matrix has eigenvalues 
(—3,-3.13,0.13). The first corresponds to a flow along the direction d g away from the 
AdS& fixed point as we flow into the IR, the second and third describe flow away and 
towards the fixed point along two orthogonal directions inside the plane of {dh,dj} AdS 4 - 
Around the Lifshitz fixed point, the flow matrix takes the Jordan normal form 
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\ 




-4 1 


(4.6) 


-4 / 





in for example the basis 



/ i9h + \dj \ 



±<9 



\ d. 



(4.7) 



,, \dh J 



There are two relevant directions (with logarithmic mixing between them) and one marginal 
direction around this fixed point. 

One can numerically solve the system ( |4.5| ) using the "shooting" technique. Since the 
AdS fixed point has relevant perturbations, it is hard to hit on the nose by varying the 
UV perturbation. Instead, we start with the AdS fixed point close to r = and try to hit 
the z = 2 critical point by integrating the flow out towards r = oo; this is much easier to 
do, because the UV fixed point has only relevant and marginal perturbations, and is hence 
(almost) UV attractive. One can find good flows that asymptote between the two fixed 
points in this way; an example is displayed in Figure 1 below. 





Fig. 1. Plotted are respectively g(r), h(r) and j(r) as functions of log{r). Their initial values 
are set at r = 1 to be respectively v/f, 0.00728 and 0.01, i.e. small deviations away from the 
Ad Si fixed point along the unique irrelevant direction. As the figure shows, all three functions 
asymptote to 1 at large r. Numerical integration of the system ( |4.5| ) suggests that the marginal 
direction at the Lifshitz fixed point becomes relevant at the nonlinear level. 



5. Future directions 

There are many directions for further research. They range from straightforward 
extensions of this work to somewhat speculative but intriguing possibilities: 
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• It would be interesting to compute more complicated observables in these spacetimes, 
including higher-point correlation functions and Wilson loops. 

• By adding bulk "probe" fields and varying the parameters in their Lagrangian, it should 
be possible to find interesting quantum phase transitions in this system. 

• The embedding into a full string theory construction of our solutions is left as work for 
the future. We anticipate that this should be possible using standard techniques of flux 
compactification. 

• A better understanding of the global and causal structure of the spacetime ([O]) would 
also be desirable. In particular, it is necessary to understand the physics behind the 
peculiar behavior at r = in more detail. This kind of behavior has been observed before 



in string theoretic brane systems ||20|| . It is possible that in this case, the behavior is 
related to the strong infrared singularities which one might expect in interacting Lifshitz- 
like theories.S 

• Our solutions have correlation functions that do not exhibit ultra- locality at zero tem- 
perature. However, in the work of 0, it was found that the Lifshitz fixed point has 
ultra-local correlators at finite temperature T, but not at vanishing T. (Ultra-locality was 
also observed in the SU{2) gauge theories discussed by [0). It is therefore interesting 
to construct finite T analogues of our spacetime (which presumably amounts to finding 
black hole solutions), and see whether the correlators exhibit some of the same behavior 
discussed in |J. 

• Finally, in any application of gauge/gravity duality to a strongly coupled system, it is 
important to understand where the large iV degrees of freedom in the field theory are 
supposed to arise. In many field theories of interest to condensed matter theorists, there 



is of course no large iV and no 't Hooft-like expansion. However, the results of [14] suggest 
that novel 3D gauge theories with SU(N) gauge group may well have Lifshitz-like fixed 
points (though they focused on the case N = 2). It seems worthwhile to try and generalize 
their results to the planar limit of large iV theories; this would potentially give a direct 
point of contact between a class of novel field theories, and gravity in spacetimes of the 
sort we studied here. 



We thank S. Dubovsky for discussions of the peculiar infrared behavior of these theories. 
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